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It has been shown by many independent studies that the cold dark matter scenario produces 
singular galactic dark halos, in strong contrast with observations. Possible remedies are that either 
the dark matter is warm so that it has significant thermal motion or that the dark matter has 
strong self interactions. We have combined these ideas to calculate the linear mass power spectrum 
and the spectrum of cosmic microwave background (CMB) fluctuations for self-interacting warm 
dark matter. Our results indicate that such models have more power on small scales than is the 
case for the standard warm dark matter model, with a CMB fluctuation spectrum which is nearly 
indistinguishable from standard cold dark matter. This enhanced small-scale power may provide 
better agreement with the observations than does standard warm dark matter. 
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I. INTRODUCTION 



Dark matter is a necessary ingredient in the standard Big Bang model of the universe. Its presence has an impact 
from subgalactic dynamics to the global evolution of the universe. However, the nature of the dark matter remains 
unknown. So far, the cold dark matter model has been very successful in explaining how structure forms 0^]. In this 
(f") ■ model the dark matter consists of weakly interacting massive particles (WIMPs) which are extremely non-relativistic 
I when structure formation begins. Because they are so massive they do not free stream and perturbations on small 
scales are preserved. In the 1980s it was realised that CDM produces too much small-scale structure, and that 
some modification of the model is needed. Several possibilities exist: there could be a large component of hot dark 
, matter damping small scale fluctuations or there could be a non-zero cosmological constant. Recent data from type 
£1h la supernovae indeed suggest that the energy density of the universe is dominated by a cosmological constant ||. 
Thus, the problem with CDM is at first sight remedied. However, in the past few years very high resolution N-body 
' simulations of structure formation have shown that any type of CDM model produces far too much substructure on 
galactic scales, compared with observations. The halo of a galaxy like our own should contain of the order 1000 distinct 
subhaloes, a factor of ten more than is found by observations ||||. Another, related problem is that galaxies are 
predicted to have singular cores. Navarro, Frenk and White || found that N-body simulations predicted a universal 
core profile of halos where p oc r _1 . Later simulations with higher resolution find an even steeper profile [f^-|To|]. At 
the same time galactic rotation curves indicate dark matter halos with finite cores, i.e. constant core density [ pd| . 
This problem is very severe and is consistently found in all simulations. 

If the details of star formation and feedback do not solve the problem, then physics at a more fundamental level 
possibly could. One option is that the primordial power spectrum has a sharp drop at subgalactic scales so that 
substructure is prevented from forming Another option along this line is that the dark matter is not cold, but 
warm Jb|[l4|]. In this model the dark matter particle mass should be around 1 keV so that the dark matter has 
significant thermal motion and perturbations on small scales are erased. However, the cut-off scale needed for the 
correct core radius of halos to be produced is so large that it is difficult to form the correct number of dwarf galaxies 

A radically different explanation was suggested by Spergel and Steinhardt |16| , namely that the dark matter could 
be cold, but have significant self-interactions. If the mean free path of the dark matter particles is of the order the 
size of the collapsing system, then the core singularity would form much more slowly, while the outer parts of the 
halo would remain unchanged. Recently, a large number of papers have appeared which investigate this possibility 
numerically Jij] -|22|]. The conclusion is that if the interactions are very strong, the model does not fit observations 
|L8| |l],^3). The halos become completely spherical apart from a small rotational deformation, and a singular core 
develops. However, it seems that models where the dark matter mean free path is similar to the system size produce 
halos closely resembling the observed ones |lj|,[l9|]. It has also been suggested that the self-interacting matter could 
be in the form of a scalar field VM. 
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That dark matter could have self-interactions is an old idea. It was originally suggested by Raffelt and Silk |25[] that 
HDM neutrinos could have strong self interactions. In this way free streaming would be suppressed and fluctuations 
only washed out via diffusion. The scenario was elaborated on by Atrio-Barandela and Davidson p6| who did a 
numerical study of this model. The possibility of number changing self interactions has also been considered [^7|-^9|. 

In the present paper we wish to explore the possibility that dark matter has both significant thermal motion and 
self-interactions. The self-interactions are assumed to consist only of two-particle scattering. In general, the inclusion 
of self interactions leads to less small scale suppression of perturbations because the small scale cut-off in power is 
given by the Jeans scale which is smaller than the free-streaming scale. We find that self-interacting hot dark matter, 
as suggested by Refs. |2j||2(|, is clearly ruled out because it produces far too little small-scale structure. However, self 
interacting warm dark matter may be a viable possibility. Strong self interactions push the power spectrum towards 
smaller scales by roughly a factor of 1.6, which may make it consistent with observations. 



II. THE BOLTZMANN EQUATION 



The evolution of any given particle species can be described via the Boltzmann equation. Our notation is identical 
to that of Ma and Bertschinger (MB) |50). We shall work in synchronous gauge because the numerical routine for 
calculating matter and CMB power spectra, CMBFAST |31{ ], is written in this gauge. As the time variable we use 
conformal time, defined as dr = dt/a(t), where a(t) is the scale factor. Also, as the momentum variable we shall use 
the comoving momentum qj = apj . We further parametrize qj as qj = qrij , where q is the magnitude of the comoving 
momentum and rij is a unit 3-vector specifying direction. 

The Boltzmann equation can generically be written as 

W\ = % = c[f], (i) 

where L[f] is the Liouville operator. The collision operator on the right-hand side describes any possible collisional 
interactions. 

We then write the distribution function as 



f{x\q,n 3 ,T) = / (g)[l + *(:r\g,n i) T)], 



(2) 



where fo(q) is the unperturbed distribution function. For a standard fermion which decouples while relativistic, this 
distribution function is simply 



f Q (q) = [eM<l/To) + l]-\ 



(3) 



where To is the present-day temperature of the species. For a self-interacting species in scattering equilibrium the 
distribution is instead 



/„(</) = [exp((e-M)/aT) + l]- 1 , 



(4) 



where e = \/ q 2 + a 2 m 2 and fi is a chemical potential. This distribution is in general different from the one for 
collisionless particles, so that one might worry that a detailed calculation of fo(q, t) is needed. However, the relevant 
quantity to look at for our purpose is the entropy per particle, s/n, which is conserved for both interacting and non- 
interacting species (note that this would not hold in a model with number-changing self interactions p7|-p9|). This 
means that for instance (p/T 7 ) = constant. Thus we do not need to worry about how the unperturbed distribution 
is changed by self-interactions. In practise we just assume that the distribution function is equal to what it would be 
for a collisionless species. 

In synchronous gauge the Boltzmann equation can be written as an evolution equation for \E' in fc-space [p0[ 
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L if\ = -H- + + -j; 

or e a In g 



h + 6r] 2 

n — m 



fo 



C[f], 



(5) 



where \i=.m?kj. h and rj are the metric perturbations, defined from the perturbed space-time metric in synchronous 
gauge Q 



ds 2 = a 2 (r)[-dT 2 + (Sij + h ij )dx i dx j ], 



(6) 
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h l2 = J d 3 ke tls (kikHk, t) + (kk ~ \kj)^{k, r)j . (7) 
Collisionless Boltzmann equation — At first we assume that j^C[f] = 0. We then expand the perturbation as 

oo 

* = £(-i)'(2l + l)^(/i). (8) 

1=0 

One can then write the collisionless Boltzmann equation as a moment hierarchy for the ^/ by performing the angular 
integration of L[f] 

*o = -k^x + 7 h—^ 9 
e 6 a In q 

*i = fc|-(* -2* 2 ) (10) 

* 1 _*A<„ 1 _ w , ) _Q. 4 + §*)^ (ii) 

= fc ( 2 ^ + l) f ( ^- 1 ~ (/ + ' 1 ~ 3 (12) 

It should be noted here that the first two hierarchy equations are directly related to the energy-momentum conservation 
equation. This can be seen in the following way. Let us define the density and pressure perturbations of the dark 
matter fluid as |3fj] 

S = Sp/p (13) 
e^ikjST^/ip + P) (14) 

<T = -{hk-\Zii){T ij -<5 ij 'r*/3). (15) 

Then energy and momentum conservation implies that j3Cj 

i =-< i+ ")( e+ i)-4(^-) i < w » 

9 = *(l-3a,)»- — 9 -*V (17) 
a 1 + w 1 + w 

By integrating Eq. @ over q 2 edq, one gets Eq. ( plf ) and by integrating Eq. ( |l(i| ) equation over c^dg one retrieves 
Eq. ©. 

Collisional Boltzmann equation — We now introduce interactions by lifting the restriction that 4jC[/] = 0. Ideally, 
one should calculate the collision integrals in detail for some explicit interaction. However, we shall instead use the 
cruder, but more model independent relaxation time approximation. Here, the right hand side of the Boltzmann 
equation is in general written as |32| 

1 * 

T C[f] = — , (18) 

JO T 

where r is the mean time between collisions. However, in this simple approximation we run the risk of not obeying 
the basic conservation laws. The collision term in Eq. (^j) is J dflj^C[f] and the one in Eq. ( |l0| ) is j d£lfij^C[f]. 
Integrating these two terms over momentum space one gets the collision terms in Eqs. (|l6)-|l7j) to be 

C[f]dn q 2 dqe (l')l 

and 

' C[f]dnq 2 dqpq = k l I C[f]dnq 2 dqq t (20) 
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respectively. However, any integral of the form 



C[f]dn q 2 dqA, 



(21) 



where A 6 (/, e, g,) is automatically zero because A is a collisional invariant (however, conservation of particle number 
(I) only applies to 2 <-> 2 scatterings). Thus, both the above integrals are zero, and the right hand side of the I = 
and 1 terms should be zero, reflecting that energy and momentum is conserved in each interaction. Apart from these 
two terms we put 



1 * ( 

yC[/]i>2 = . 

JO T 



so that the full Boltzmann hierarchy, including interactions, is 



, q T 1 ; d In /o 
-k^x + -h— j-^ 
e b fling 

fc|-(* -2* 2 ) 



fcf (2*! -3* 3 )- I ^h+^f, 



5e 
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(2Z + l)e 



(Z*i_i - (Z + l)*,+i) 



dln/o * 2 
din g r 

T 



Z > 3 



(22) 

(23) 
(24) 
(25) 
(26) 



In Appendix A we discuss how the above set of equations relates to the equations used in other studies of self-interacting 
dark matter. 

Relaxation time — We now need an expression for the collision time r. In general we can write 



t 1 = 7T,(tr|w|). 



(27) 



For relativistic particles scattering via exchange of a massive vector boson (mx 3> T, m where mx is the vector boson 
mass and m is the mass of the dark matter particle) we have 



whereas for non-relativistic particles it is 



(a\v\) cx (T/to) 2 , 
(a\v\) oc (T/m) 1 / 2 . 



As an interpolation we use 



<aM> 



-0-0 



1/2 



(28) 
(29) 

(30) 



III. NUMERICAL RESULTS 



Using the above equations we have calculated matter and CMB power spectra for two different dark matter models: 
HDM (m = 10 eV) and warm dark matter (m = 1 keV) over a range of scattering cross sections. In practice we have 
incorporated the equations into the CMBFAST code developed by Seljak and Zaldarriaga pljj . All the models were 
done assuming that fix = 0.95 and Qb = 0.05, Hq = 50kms _1 Mpc -1 . The conclusions are unchanged if a ACDM 
model is assumed, since our purpose here is only to show how self-interactions change the power spectra. Fig. 1 shows 
the matter power spectrum in terms of the quantity 

A 2 (fc) ^ (31) 

for our two different cases. In both cases, the power spectrum cut-off is pushed towards higher k if self-interaction is 
assumed. The HDM (m = 10 eV) results are in agreement with the results of Atrio-Barandela and Davidson (26), for 
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k smaller than the cut-off scale. At small scales, their results are somewhat different from ours, probably because of 
an erroneous term in their perturbation equations (as explained in the appendix). 

For our choice of particle masses, the dividing line between the non-interacting and strongly interacting regimes is 
roughly at 

a ^ l(T 36 cm 2 . (32) 




FIG. 1. Matter power spectra for two different dark matter particle masses. The dashed line is for no self-interaction, the 
triple-dot-dashed is for <jo = 1.2 x 10~ 38 cm 2 , the dot-dashed for <to = 8 x 10 _36 cm 2 , and the dotted is assuming complete 
pressure equilibrium. For reference we have plotted the spectrum for standard cold dark matter (full line). 

Note that this is much lower than the cross section which is needed to explain structure on galactic scales in the 
self- interacting cold dark matter model. In that case, the dividing line is closer to 10 -23 cm 2 . For the case where the 
dark matter is hot, self-interactions are not able to improve the agreement with observations significantly because 
the power spectrum cut-off is still at much too large a scale. As discussed in Ref. |jq] , warm dark matter provides 
a good fit to observations of dwarf galaxies if the power spectrum cut-off is at roughly 2h 50 Mpc~ 1 , corresponding 
to a mass of 1 keV. However, explaining the core structure of dark matter halos requires that m < 300 eV [fl5| , so 
that even though the uncertainties involved in determining the best cut-off scale are as large as a factor two [ fL5l , the 
collisionless warm dark matter model is inconsistent with observations. Our results indicate that it might be possible 
to lower the warm dark matter particle mass to this smaller value and compensate by making the warm dark matter 
self-interacting, which decreases the cut-off length scale by about a factor of 1.6 compared to the non-self- interacting 
case. Numerically we find that the k where A 2 (fc) takes its maximum value is well approximated by 
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A 2 (fc)max _ f 1 - 1 (tEv) 3/4 mPC 1 collisionless, 

\ 1.7 ( ik^v ) 3 ^ 4 Mpc -1 strongly self-interacting. 

For the collisionless case this corresponds to the free-streaming scale, whereas in the strongly interacting case it 
corresponds to the Jeans scale for a given particle mass. From this result we conclude that self-interacting warm dark 
matter is marginally consistent with the present observational constraints. 

For the CMB, the fluctuations are usually expressed in terms of the Ci coefficients, C; = (|ai m | 2 ), where the a; m 
coefficients are determined in terms of the real angular temperature fluctuations as T(6, 4>) = ^2 im ai m Yi m (9, <fi). Fig. 2 
shows the CMB spectra for the same two particle masses. If the dark matter is hot, the CMB spectrum is changed 
relative to cold dark matter, because the DM particles are not completely non-relativistic at recombination. This 
gives rise to what is called the early integrated Sachs- Wolfe (ISW) effect. Self-interactions have very little impact 
because they only affect scales within the dark matter sound horizon at recombination. Even for a dark matter mass 
of 10 eV, this is at too small a scale to have a significant impact. For a dark matter particle mass of 1 keV, the effects 
are completely negligible. Our results for non-self-interacting warm dark matter agree with those of Burns |35| ]; we 
have extended his results to demonstrate that the addition of self-interactions to the warm dark matter model also 
produces a negligible difference from standard CDM. 
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FIG. 2. CMB power spectra for the same models as in Fig. 1. The curve labels are also identical to those in Fig. 1. 



IV. DISCUSSION 



We have performed a quantitative calculation of the linear behaviour of warm dark matter models with possible 
self interactions. As expected, power on small scales is generally increased in self-interacting models because free 



6 



streaming is suppressed. In collisionless models, power is suppressed on the free streaming scale, whereas in strongly 
self-interacting models the cut-off is at the Jeans scale. This increase in the amplitude of the fluctuations on small 
scales has the effect of pushing the cut-off in the power spectrum down to smaller scales by approximately a factor of 
1.6. This may allow warm dark matter to better fit the dwarf galaxy observations for masses which are small enough 
to explain the core structure of dark matter halos, a result which could make warm dark matter a more viable dark 
matter candidate. 

Our CMB results indicate that, like standard warm dark matter, self-interacting warm dark matter is indistinguish- 
able from standard cold dark matter in terms of the CMB fluctuation spectrum. Thus, it is one of the few variants 
on the standard model which will not be probed by future CMB experiments. Any constraints on this model must 
therefore come from large-scale and galactic structure considerations. For instance, analysis of high-z structure like 
damped Ly-a systems might lead to interesting constraints. 

Note that the cross section for scattering of dark matter particles would have to be of the order 10 -36 cm 2 in order 
to change the matter power spectrum significantly. This is orders of magnitude more than the cross sections typical 
in weak interactions, and at present there are no obvious candidates for such dark matter particles. However, it could 
well be that warm dark matter with relatively strong self-interactions could be in a mirror sector, in which case there 
are no real restrictions 1341 . 
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APPENDIX A: THE BOLTZMANN EQUATION IN DIFFERENT ASYMPTOTIC LIMITS 



1. Large scattering cross sections 



In the limit of very large scattering cross sections, the dark matter is kept in pressure equilibrium until the present. 
This is the type of evolution assumed in Refs. |27] In this case the evolution equations read 

*o = -***i + s^ (AD 
e 6 ding 

*i = fcftfo (A2) 
* ; > 2 = 0. (A3) 
By performing the appropriate momentum integrations this yields 

5=-(l + cj)(d+^] 3- (A4) 




= ~(l-3w)0 9+ — '—k 2 S (A5) 

a 1 + u 1 + lo 

This equation is equivalent to Eqs. (13-14) in Ref. [^9| (when their T = II = 0), which are written in gauge invariant 
form. 



2. Large k limit 



At very small scales one may as a first approximation neglect the metric perturbations. The Boltzmann hierarchy 
can be truncated by neglecting terms higher than second order (including a) , similar to how the Enskog expansion is 
performed |32] ] . Then the hierarchy equations when integrated over momentum yield 

6 = /j 2 (<5/4-46»t/15). (A6) 
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It is interesting to compare our set of equations with Eqs. (25-26) of Atrio-Barandela and Davidson (AD) (g6[. They 
are almost identical, except for the term proportional to H in their equation. For relativistic particles this term should 
be zero, as it is in the above equation. 

The term 40t/15 can be interpreted as a shear viscosity term, which can in general be written as r]9/p prji . Here 
r\ is the viscosity of the fluid. Using this parametrization we find that 



For a relativistic gas with Boltzmann statistics, p = 3Tn, so that 



V 



—Tut. 

5 



(A7) 



(A8) 



This expression for the fluid viscosity agrees with what is found in Ref. |26| (their Eq. (33)). From Eq. A£, one can 
see that the perturbations oscillate and are damped at the rate 

2 ' (A9) 



r = — rk 2 . 
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